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We report on high-resolution interferometer observations of the ^^CO(J=1-0) emission in the central 5 
kpc region of the Seyfert 2 galaxy NGC 4501. The observations were made using the Nobeyama Millimeter 
Array during a long-term CO line survey of Virgo spirals. The major features are: (1) a nuclear concen- 
tration with a radius of r ~ 5" (390 pc), which is resolved into double peaks, and (2) spiral arms which 

^vq ' extend out from the nuclear region. The nuclear component has a mass of 1.3 x 10® M©, which corresponds 

to ~3.5% of the dynamical mass, and shows a slight non-circular motion. The double peaks are separated 
by ~ 4."7 (370 pc), and located on the root of optical spiral arms in a HST image. The gas arms are 
^ I associated with the spiral dust lanes, and are linked to the central double peaks. The non-circular motions 

r^ . along the molecular arms indicate the fact that the gas is driven by the density wave, rather than the 

stochastic processes. We calculated the gas cloud orbits in a stellar spiral potential, and explained the 
observed CO spiral arms and non-circular motions. We suggest that the central gas condensation arises 

1^ ■ from spiral-driven gas transfer. We estimated and compared the effect of two possible mechanisms of 

t:;j- \ angular-momentum transfer: galactic shock, and gravitational torques. We discuss that the galactic shock 

^^ ■ is dominant. 

^H \ Key words: galaxies: individual (NGC 4501) — galaxies: ISM — galaxies: kinematics and dynamics 

O . — galaxies: Seyfert — galaxies: spiral 

o 



1. Introduction in full flourish several decades ago (e.g. Fujimoto 1968b; 

^ , Roberts 1969), its application for observational data and 

Gas dynamics in disk galaxies is closely related to discussions within the context of radial mass transfer has 

the stellar potential. Gas dynamics in barred galax- not been thoroughly pursued compared to that in barred 

ies has been well studied in particular, and many ob- galaxies. 



?H \ served features have been successfully explained. Bars We observed a non-barred CO-luminous galaxy, NGC 

^ are thought to be the most viable candidates for fuel- 4501, during the Virgo high-resolution CO survey (ViCS: 

ing the central activity, transferring mass from large to Sofuc et al. 2003), using the Nobeyama Millimeter Array 

small scales (e.g. Simkin et al. 1980; Shlosman et al. 1989; with the highest resolution ever achieved for this galaxy. 

Athanassoula 1992). Recent statistics report circumstan- The prominent molecular arm structure, non-circular mo- 

tial evidence of radial gas inflow as a result of bar-driven tion, and centrally-condensed double peaks motivated our 

gas dynamics. Metallicity gradients in barred spirals are work on the gas dynamics in the spiral potential and pos- 

systematically flatter than that of unbarred ones, and they siblc inflow mechanism to the center, 

become shallower when the bar length and the bar ellip- NGC 4501 is one of some non-barred galaxies that have 

ticity increase (Martin, Roy 1994). The degree of gas a relatively high degree of central gas condensation, simi- 

condensation to the central kiloparsec is higher in barred lar to those of barred galaxies (Sakamoto et al. 1999b). It 

galaxies than in unbarred ones (Sakamoto et al. 1999b). is an SAb galaxy (de Vaucouleurs et al. 1991, hereafter 

However, the central gas condensation is also prevalent RC3) at a distance of 16.1 Mpc (Ferrarese et al. 1996; 

in non-barred galaxies (Sakamoto et al. 1999a), and some 1" corresponds to 78 pc) located about 2.°0 (0.56 Mpc) 

studies claim that the level of the central activity does northeast of the center of the Virgo cluster, M 87. It 

not correlate with the presence of large-scale bars (Ho et hosts a weak Seyfert 2 nucleus (Ho et al. 1997a). It is 

al. 1997b; Mulchaey, Regan 1997). Their results suggest slightly Hi deficient (Chamaraux et al. 1980), and has a 

that gas inflow may occur both in barred and non-barred high CO luminosity and a high IRAS flux density (Stark 

spiral galaxies. Although a theoretical study of the gas et al. 1986). While it is not classified as a grand-design 

dynamics in spiral galaxies with tightly wound arms was spiral galaxy in i3-band image (arm class 9; Elmcgrcen, 
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Table 1. Parameters of NGC 4501. 



Table 2. Observation Parameters 



Parameter 



Value 



Source* 



Hubble type 
Nuclear activity 
Adopted distance (Mpc) 
P.A.(isophotal)(°) 
Inclination(°) 

Vsys (kms^i) 
^t°otai (mag) 
5co(45")(Kkms-i) 
Linear scale (pcarcsec^^) 



SA(rs)b 1 

Type 2 Seyfert 2 

16.1 3 

140 4 

58 4 

6.'9 X 37 1 

2263 4 

9.86 1 

12.0±4 4 

78.1 



* (1) de Vaucouleurs et al. (1991) (RC3); (2) Ho et al. (1997a); (3) 
Fcrrarcse et al. (1996); (4) Kenney and Young (1988) 



Elmegreen 1987), it shows continuous two-armed spiral 
in the iT-band within r ^ 3kpc (Elmegreen et al. 1999, see 
their figure 1). Moreover, CaroUo et al. (1998) noted the 
spiral-like dust lanes down to the nucleus in the F606W 
(wide y-band) image of HST WFPC2. The basic proper- 
ties of NGC 4501 are hsted in table 1. 

In this paper, we show that the gas dynamics in the 
central region of this galaxy is governed by a stellar spiral 
potential, and suggest the possibility that the observed 
central features can be the outcome of radial inflow due 
to the spiral arms. High-resolution ^^CO(J=1-0) observa- 
tions of NGC 4501 are presented in section 2. The results 
from the observations are presented in section 3. The CO 
distribution and the velocity field are modeled based on 
observations in section 4. Possible mechanisms of spiral- 
induced radial gas inflow are discussed in section 5. We 
summarize our conclusions in section 6. A detailed de- 
scription of the damped-orbit model in the spiral potential 
is shown in appendix 1 . An estimation of the change in the 
angular momentum with shocks and torques is presented 
in appendix 2. 

2. Observations and Reduction 

2.1. CO(J=1~0) Observations with NMA 

We carried out observations of NGC 4501 in the 
^^CO(J=1-0) line using the Nobeyama Millimeter Array 
(NMA) at the Nobeyama Radio Observatory (NRO) dur- 
ing the course of a long-term project, the Virgo high- 
resolution CO Survey (ViCS; Sofue et al 2003). The obser- 
vations were made from 2000 December to 2002 March for 
a single pointing center at (aj2ooo, '5j20oo) = (12''31'°59?14, 
+14°25'12."9). The NMA consisted of six antennas, each 
had a diameter of 10 m, providing a FWHP of about 
65" at 115 GHz. We used the spectro-correlator Ultra- 
Wide-Band Correlator (UWBC: Okumura et al. 2000), in 
a mode of 256 channels covering 512 MHz. One channel 
(2 MHz) corresponded to 5.2kms~^ at the observing fre- 
quency. Quasar 3C 273 was observed every 20 minutes to 
correct for any instrumental gain variations and the band- 
pass response. Continuum observations of Mars were used 



Observed center frequency (GHz) 


114.407 


Array configurations 


AB, C, and D 


Observing field center 




(Sakamoto et al. 1999a): 




a(J2000) 


12i>31™59?14 


J(J2000) 


-hl4°25'12."9 


Frequency channels 


256 


Total bandwidth (MHz) 


512 


Velocity coverage (kms~^) 


1342 


Velocity resolution (kms~^) 


5.24 


Amplitude and phase calibrator 


3C273 


Primary beam (") 


65 



to calibrate the flux scale. We calibrated the absolute flux 
scale in each sideband in each observation: once in 2000, 
three times in 2001, and once in 2002. We compared the 
noises in two sidebands and used the lower sideband (102 
GHz) data in all flux calibrations. The flux density of 3C 
273 was roughly constant at 9.6 Jy, within ± ~ 20%. No 
significant flux variation was observed in all of our runs. 
The chopper-wheel method was used to correct for any 
atmospheric transmission losses and elevation-dependent 
gain variations. We used AB, C, and D array configura- 
tions of the NMA to achieve a resolution of less than 2" . 
The visibility data cover projected baselines from 10 to 
341 m. Partially shadowed observation data were deleted, 
and the minimum baseline length projected onto the sky 
was the diameter of each antenna. This resulted in the 
central hole in the u-v plane, which restricts the largest 
detectable size of objects to about 54". The observation 
parameters are listed in table 2. 

2.2. Reduction 

The raw visibility data were calibrated with the 
NR0/UVPR0C2 software (Tsutsumi et al. 1997) and 
mapped with the NRAO/AIPS package. We applied the 
CLEAN procedure with two ways of weighting (natural 
and pure uniform), and obtained three-dimensional data 
cubes (R.A., Decl, V^sr) in low (5") and high (2") spa- 
tial resolutions. The synthesized beam parameters, veloc- 
ity resolutions and rms noises for these maps are listed 
in table 3. The fraction of the total flux recovered by 
our NMA data cube was estimated by comparing with 
single-dish observations using the FCRAO 14 m telescopes 
(Kenney and Young 1988). For this comparison, the data 
cube was corrected for the primary beam attenuation of 
the NMA antennas, convolved with the single-dish beam 
(FCRAO, 45"), and sampled at the pointing centers of the 
single-dish observations to obtain a value comparable to 
the single-dish flux. The NMA cube was found to recover 
about 78% of the single-dish flux in the central 45". Note, 
however, the primary beam correction was not applied to 
the maps shown in this paper. 

Figures 1 and 2 display low- and high-resolution inte- 
grated intensity maps and velocity fields in the central 
region of NGC 4501 (90" ~ 7.0 kpc square for figure 1, 



No. 



Gas Dynamics in NGC 4501 

Table 3. Parameters of Maps.* 





Weighting 


Beam 


Velocity 
Resolution 

(kms-i) 


rms noise 

a 

(mjybeam^^) 


Tbfor 


Resolution 


FWHM P.A. 

(") n 


1 Jybeam"^ 

(K) 


Low 
High 


Natural 
Pure uniform 


5.6 X 3.7 160 
1.8x1.7 138 


10.5 
41.9 


14 
16 


4.5 
30.5 



* The map centers are set at derived dynamical center: 

(aj2000, <5j20oo)=(12''31™59.n2, +14°25'13."3) and derived systemic velocity 2261 kms"^ 



14" ^ 1.1 kpc for figure 2). The crosses in figure 1 and 
figure 2 indicate the major and minor axes of NGC 4501 
(P.A. =140°). The ellipse in figure 1 indicate a projected 
circle of diameter 75" (5.9 kpc), inclined with the position 
angle and inclination (58°) of this galaxy. The size of the 
crosses in high-resolution maps is 5" (390 pc). The cen- 
ter of all crosses and ellipses in figures 1 and 2 are set at 
the dynamical center derived from the velocity field (see 
below) . 

Figure 3 shows channel maps with an interval of 4 
MHz (10.4 kms"-'^) in the same region. Significant emis- 
sion (> 3ct, 1(7 = 14mJybeam^^) was detected in 46 
adjacent channels within the velocity range of V lsr. ~ 
2044-2513kms-i (A^ = 269kms-i). In order to check 
the contribution of the continuum flux, we examined the 
lower sideband (102 GHz) data that can be regarded as 
being emission-free. We summed up our data over the 
full bandwidth (512 MHz); no significant continuum flux 
(> 3(7, 1(7 = 1.3mJybeam^^) was detected. 

We obtained kinematical parameters from the low- 
resolution (5") velocity field using the AIPS/GAL pack- 
age. The dynamical center (table 2) was determined using 
the Brandt rotation curve model (Brandt, Scheer 1965). 
For fitting, we used the central r < 5" region where the 
velocity field can be regarded as being nearly symmetric 
against the dynamical center. The initial guess of the pa- 
rameters affected the result within r < 1" . The centers of 
CO maps and PV diagram are set at this derived dynam- 
ical center: (ajaooo, '^J20oo) = (12''31"^59?12,+14°25'13"3) 
and the systemic velocity T^ys — 2261 km s~^. A GAL fit 
also gave the inclination and position angle of the inner 
disk. However, these quantities are dependent on the non- 
circular rotation, which are not negligible in the central 
regions. Therefore, we adopted the inclination (58°) and 
position angle (140°) determined from optical images of 
the entire stellar disk. 

Figure 4 shows the radial profile of the CO-line inten- 
sities in units of Jy kms~^ arcsec"^, which is corrected 
for the inclination. At first, low-resolution integrated- 
intcnsity map was corrected for primary-beam attenu- 
ation. We then applied the AIPS task IRING to az- 
imuthally average the integrated intensities in each an- 
nulus with corrections for the inclination. We applied the 
central positions derived above and the inclination and 
position angle provided from optical observations (table 
1). Figure 5 displays a position- velocity diagram (PV di- 
agram) along the optically defined major axis (P. A. =140°) 
with a slit width of 5" and 60", respectively. The veloci- 



ties have not been corrected for the inclination. 

2.3. Supplementary Data 

We obtained optical photographs of NGC 4501 from 
archives for a comparison. We overlaid low-resolution CO 
map contours on a Digitized Sky Survey _B-band image 
in order to see the global pattern (figure 6). We also 
overlaid high-resolution CO map contours on an unsharp- 
masked image made from the HST WFPC2 observations 
at the F606W (wide F-band) filter (Carollo et al. 1998), in 
order to see the innermost r < 500 pc region with a higher 
resolution (figure 7). 

3. Results 

The obtained CO map (figure 1) shows spiral arms and 
a strong central condensation, the central CO condensa- 
tion is resolved into double peaks in the higher resolution 
map (figure 2). We describe these components in the fol- 
lowing subsections. 

3.1. Main Disk 

3.1.1. Distribution in the main disk 

The entire CO disk of NGC 4501 has an exponential 
profile of scale length 0.'7, and it continuously extends to 
r ^ 3f8 (Kenney and Young 1988). We now discuss the 
central r ^40"(3.1kpc) region of the CO main disk (figure 
1, Left). The radial profile is nearly exponential (figure 
4), and the emission enhancements due to the molecular 
spiral arms are distinct (subsection 3.2). 

3.1.2. Kinematics of the main disk 

The velocity field (figure 1, Right) shows rotation of 
the main disk with the north-west side approaching and 
the south-east side receding. It represents non-circular 
motions superposed on regular rotation of the disk. The 
non-circular motions appear as a distortion of the isove- 
locity contours. The isovelocity contours around the spi- 
ral arms of r ^ 2.9 kpc (the ellipse in figure 1) indicate 
that the near (northeastern) sides of the arms are reced- 
ing, and the far (southwestern) sides are approaching, as 
compared with pure circular motions. The non-circular 
motions also appear in r ^ 10" as a Z-shaped pattern of 
the isovelocity contours. They can be attributed to spiral 
arms (subsection 3.2). The channel maps show the rota- 
tion of the disk and non-circular motion along the spiral 
arms as well (figure 3). The map centers were set at the 
derived dynamical center (section 2). Figure 9 shows the 
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Fig. 1. Low-resolution (5."6 X 3."7) CO(l-O) intensity map and velocity field of NGC 4501. The synthesized beam is shown in the 
lower left-hand corner. The primary-beam attenuation is not corrected. The clip level was Aa {la = 14mjybeam~^ in each channel) 
to make these moment maps. The circle represents the FWHP of NMA (65"). The crosses indicate the orientations of the major and 
minor axes of NGC 4501 (P.A.=140°). The center of maps, crosses and ellipses are set at the dynamical center (table 2). The ellipse 
indicates a projected circle of a diameter of 75"(5.9kpc) inclined with the position angle (140°) and the inclination (58°). Left: 
Intensity map in the central 90" (~ 7.0 kpc) squared region. The contours are drawn at 5, 10, 20, 30, 40, 60, 80% of the peak value 
46Jybeam~^ kms~^. The rms noise is Icr = 1.0 Jybeam"^ kms~^. 1 Jybeam"^ corresponds to 4.51K. Right: Velocity field. The 
contours are drawn every 50kms^^ in the range of 2050— 2550 km s"'^. The white line indicates the systemic velocity, 2261 kms~^. 
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Fig. 2. High-resolution (l."8) CO(l— 0) intensity maps and velocity field of NGC 4501. The synthesized beam is shown in the 
lower left-hand corner. The primary-beam attenuation is not corrected. The clip level was 4cr (Icr = 16mjybeam~^) to make 
these maps. The crosses indicate the dynamical center and the orientation of the galactic major and minor axes. Their size are 
5"(390pc). Left: Intensity map in the central 14"(~ 1.1 kpc). The contours are drawn at 30, 40, 50, 60, 70, 80, 90% of the peak 
value 9.6 Jybeam"^ kms~^. The rms noise is Icr = 2.1 Jybeam"^ kms~^. IJybeam"^ corresponds to 30. 5K. Right: Intensity 
map representing velocity field, in which approaching {V < Ksys: solid lines) and receding (V > Vsys: broken lines) components are 
distinguished. The contours are drawn at 30, 45, 60, 75, 90% of the peak value. 
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Fig. 3. Channel maps in the CO(l-O) emission. The elhpse is the same as that of figure 1. The channels have an interval of 
10.4kms~^. Their central velocities (Vlsr in kms~^) are labeled at the upper left corners. The contour levels are -3, 3, 6, 9, 12, 15 
a, where la = 14mjybeam~^. The negative contours are dotted. No primary-beam correction was applied. The velocity distortion 
along the spiral arms are apparent. 
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inclination-corrected rotation curve. The rotation curve 
was derived from the position-velocity diagram along the 
major axis with a slit width of 3" using an iteration 
method developed by Takamiya and Sofue (2002). This 
method determines a rotation curve so that it can repro- 
duce the observed position-velocity diagram. 

The rotation velocity shows a steep rise in the central 5" 
to w ^ 200 kms""'^, and then gradually rises to ^ 290kms~^ 
at r ^^ 35". This velocity is similar to the rotation velocity 
of the outer flat part (Nishiyama et al. 2001), but slightly 
excesses the Hi rotation velocity, 268kms~^, at 20% of 
the peak flux (Cayatte et al. 1990). 

3.1.3. Mass and surface brightness of the molecular gas 
The mass of the molecular gas M^^ is estimated from 
the total CO-hne flux, Sco^ adopting a Galactic Nn^/Ico 
conversion factor, Xqo^ a-nd the galactic distance, D, as 

2 






7.2x10^ 



D 

Mpc 



5, 



CO 



\ Jykms 



1.8 X 102"cm-2(Kkms' 



(1) 



The gas surface densities are independent of the adopted 
distance to the galaxy. If we assume the hydrogen mass 
fraction 0.707 (Dappen 2000), the total gas mass including 
other elements, such as helium, becomes Mgas = l.ilM^^. 
The surface mass density of molecular hydrogen on the 
galaxy plane is calculated from the integrated CO-line in- 
tensity, /cOi the galaxy inclination, i, and the conversion 
factor, XcO) as 



5^H2 



Mqpc- 



=3.0 X lO^cosi 



'CO 



Jykms ^ arcscc ^ 
-'^co 



1.8 X 1020cm-2(Kkms-i) 



(2) 



We adopt the conversion factor of Sco = N-n-^/Ico = 1-8 x 
1020cm-2(Kkms"i)"i (Dame et al. 2001) from observa- 
tions in the Galaxy, while Xqq could be smaller in galac- 
tic center regions (Arimoto et al. 1996). The total flux 
within r = 5kpc (32") is 5co = 3.9 x 10^ Jy kms"\ which 
corresponds to Mgas = 1.0 x lO^M© for the galaxy dis- 
tance D = 16.1 Mpc. The dynamical mass is derived from 
rotation curve as A/dyn = fw^/G, assuming circular rota- 
tion. The ratio of the dynamical mass, 4.3 x 10"'^° Mq (v = 
270kms~^), to the total molecular gas mass, Mgas/Afdyn, 
is about 2.4% within r ^ 5kpc. 

3.2. Spiral Arms 

Two molecular spiral arins extend out from the center 
to the end of our field of view in the map. The arms 
start from the northern side and the southern side of the 
nucleus, and continue to the northeast arm and the south- 
west arm, respectively. Their outer parts are less clear ow- 
ing to the primary beam attenuation. These spiral arms 
are distinct up to r ^ 50" in larger spatial coverage, but 
lower resolution map by Wong and Blitz (2002). Our re- 
sult shows these arms more clearly, owing to the higher 
resolution (~ 5"). It also shows many giant molecular 
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Fig. 4. Radial distribution of the molecular gas surface den- 
sity in NGC 4501 with a logarithmic scale. The CO in- 
tensities were corrected for primary beam attenuation and 
inclination, and averaged over annuli with 0''5 = 39.0 pc. 
Xco = 1.8 X 1020cm-2(Kkms-i)-i and Sgas = 1.41Sh2 was 
adopted to calculate the total gas surface density in unit of 
mass per parsec square. The two horizontal bars show the 
beam size along the major and minor axes (5."6x 3."7). 



cloud associations in the arms. These arms are coincident 
with the dust lanes (r ~ 2.9 kpc) in an optical image (fig- 
ure 6) . The dust lanes are continuous down to the nuclear 
region of ~ 500 pc. The line-of-sight velocities along these 
arms show deviations from that of a circular-rotating disk. 
It is about ~ 50kms^^ on the galactic minor axis. They 
appear in the position-velocity diagram as a "figure-of- 
eight" pattern (figure 5). These distortions in the veloc- 
ity field is similar to that of the Hi spiral arms of M 81 
(Visser 1980), which are typical in gaseous spiral arms 
caused by density waves. 

3.3. Central Double Peaks 

The low-resolution map shows a central condensation of 
radius ^ 5"(~ 390 pc). The total CO flux of this compo- 
nent is Sco = 48Jykms^^, which corresponds to Mgas = 
1.3 X 10^ M0. The gas-to-dynamical mass ratio within 
r = 5" (390 pc) is ~ 3.5%, since the dynamical mass de- 
rived from V ^ 200kms-i at r = 5" is Mdyn = 3.6 x IO^Mq. 
The low-resolution (5") velocity field indicates a slight 
non-circular motion. This nuclear component shows high 
velocities of ~ ±200 kms^^ in the position-velocity dia- 
gram, and is resolved into two bright peaks (figure 5) with 
distinct velocities. 

The nuclear concentration is also resolved into double 
peaks in the high resolution (^ 2") intensity map, sepa- 
rated by 4"7 (~ 370 pc) (figure 2). Each of these corre- 
sponds to the double peaks in the P-V diagram. The net 
flux within r ^ 5" is 14Jykms^^ in high resolution, while 
it is 48 Jykms"^ in low resolution, because the higher reso- 
lution map missed the extended fluxes. The double peaks 
may be on more diffused components, which are unseen 
in our image. 

The peak intensity in the high-resolution map is 
9.6 Jybeam~^kms~^ or 290 Kkms~^, which corresponds 
to Ehs = 440 Mqpc-2 and Egas = 620 Mopc'^. 
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Fig. 5. Low-resolution position-velocity diagram along the major axis of NGC 4501. Velocities are not corrected for inclination 
(58°). Left-hand side corresponds to the northwest. Contour levels are 20, 30, 40, 50, 60, 70, 80, 90% of the peak value. Left: The 
slit width is central 5". The peak value is 4.9Jybeam~l or 51 K. The central steep rise is resolved into double peaks. Right; The 
slit width is 60", including extended spiral arms. The peak value is lljybeam"^ or 110 K. The figure-of-eight pattern due to the 
spiral arms are apparent. 



4. Gas Dynamics in NGC 4501 

As we have shown above, the observed velocity field 
demonstrates distinct non-circular motions along the spi- 
ral arms. Its well-ordered pattern indicates that the den- 
sity wave, rather than the stochastic process, dominates 
in the central region of this galaxy. Although this galaxy 
seems to be multiple-armed in the outer region, there is 
some evidence that it has two arms in the central r ^ 40" . 
Elmegreen ct al. (1999) showed a continuous two-armed 
spiral in the X'-band image, and concluded that in this 
galaxy there are density waves caused by the spiral arms of 
an older stellar population. An unsharp-masked 2MASS 
ii'-band image is shown as figure 8. 

Optically unbarred galaxies sometimes unveil their hid- 
den bars in the images at longer wavelengths (e.g., Scoville 
et al. 1988; Mulchaey et al. 1997; Kohno et al. 2003), be- 
cause contributions from young stars and disturbances by 
dust lanes may make a spiral-like appearance at shorter 
wavelengths. Several studies have examined this galaxy 
for a bar, by means of elliptical isophote fitting on op- 
tical and near-infrared images (Jungwiert et al. 1997; 
Elmegreen et al. 1999; Sil'chenko et al. 1999). If a bar 
structure is dominant, the position angles, and eccentric- 
ity of the isophotal ellipses are expected to vary with the 
radius, since the bars are triaxial and their eccentricity 
shows a radial variation. Near-infrared images are prefer- 
able for this purpose because of their advantage to reflect 
a galactic potential better than optical images. The cri- 



teria for bar identification in general is: (1) the ellipticity 
increases as a function of the radius, and then decreases 
to reveal the inclination of the disk; (2) the position an- 
gle is constant over the radii where the ellipticity is rising 
(McLeod, Rieke 1995; Mulchaey et al. 1997). Jungwiert 
et al. (1997) and Elmegreen et al. (1999) independently 
conducted this analysis in the ff-band and the K'-h&nd, 
respectively, and obtained quantitatively almost the same 
results. Interior to r ~ 15", position angle remains con- 
stant at about 140°, and the ellipticity gradually rises 
from zero to ^0.4. Then, P. A. declines to reach the mini- 
mum value of ^ 120° at r '--^ 27", and gradually recovers to 
attain its original value of 140° at r ~ 50", and then keeps 
this value constantly to the outer region. The ellipticity 
also drops at r ^ 20" to 0.3, and then gradually increases 
to 0.5 at r ~ 40". (see figure 4 in Jungwiert et al. and 
figure 5 in Elmegreen et al.) 

From their quantitatively same results, Jungwiert et al. 
concluded that there is no evidence for a bar or other tri- 
axiality, while Elmegreen et al. concluded the opposite. 
Giving a certain criteria for presence of a bar, and based 
on their careful simulations of the effect of projection on 
an identification of a bar, Jungwiert et al. concluded that 
the change of ellipticity in the inner region and the con- 
stancy of the P. A., as well as the fact that its value cor- 
responds to that of the outer disk, is compatible with a 
spheroidal shape of the bulge. The slight change of the 
eccentricity and P.A. behind r ^ 18" should be attributed 
to spiral arms. 
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Fig. 6. Comparison of the low-resolution CO map (contours) and the f?-band image (gray scale, from the DSS archive). Northeastern 
CO spiral arms trace dark dust lanes. 




Fig. 7. Comparison of a high-resolution CO map (contours) and an unsharp-masked HST wide V-band image (gray scale). CO 
double peaks are located on the roots of spiral dust lanes. 
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Fig. 8. Deprojected overlay of an unsharp-maskcd A'-band 
image (gray: from the 2MASS archive). The vertical axis is 
taken to be galactic major axis. 

On the other hand, Elmcgreen et al. noted the presence 
of a small bar of 15" diameter (see their table 1), without 
any explanation about the criteria to identify the bar. It 
seems that they might judge from the growing ellipticity 
in the inner region, however, as discussed in Jungwiert 
et al.(1997), that a bar is not relevant to explain those 
characteristics in this region. 

Moreover, no large-scale bar has been found in NGC 
4501, and the HST X^-band image shows the spiral arms 
and spiral dust lanes penetrating into the nucleus (figure 
7). Based on these facts, we consider the spiral structure 
to be prominent, rather than a bar, in the gas dynamics 
of NGC 4501. 

In the following section, we calculate the gas cloud or- 
bits in a potential made by stellar spiral arms by mod- 
ifying the damped-orbit model of Wada (1994) to un- 
derstand the gas motion in NGC 4501. Based on the 
spiral-arm driven gas dynamics, we suggest mechanisms 
of the angular-momentum transfer in spiral potentials, 
which could produce the central gas condensations, and 
probably, the double peaks. 

4-.1. Damped- Orbit Model in Spiral Potential 

In order to study the gas responses in a weak spiral 
potential, we modified the damped-orbit model for a bar 
potential (Wada 1994). This model gives closed gas orbits 
in a barred potential by analytically solving the equations 
of motion, which include a damping- force term, to emulate 
the collisional nature of the gas. It has been applied to 
galaxies NGC 5005 (Sakamoto et al. 2000) and NGC 3079 
(Koda et al. 2002), and successfully explains the observed 
gas motions in these barred galaxies. Instead of a straight 
bar potential, we used a spiral potential. Modeling after 
Binney and Tremaine (1987), we assume that the stellar 
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Fig. 9. Rotation curve derived by an iteration method of 
Takamiya and Sofuc (2002) and reproduced position-velocity 
diagram by convolving the rotation curves with the observed 
intensity distribution. The velocities are corrected for incli- 
nation. 

potential is represented by an axisymmetric disk poten- 
tial, $0, and a first-order perturbation owing to the spiral 
arms, $i, as 

MR,^j)^MR) + ^iiR,^). (3) 

We use the Plummcr potential for the axisymmetric disk 
part. 

The scale length of the model is normalized by a. $i is 
the first-order perturbation by spiral arms, written as 

<Pi{R,^) = %{R)cos[m{no-^p)t+plog{R)], (5) 

which represents a logarithmic m-armed spiral with p = 
mcotO, where 6 is the pitch angle of the spiral arms. We 
set TO = 2 and 6* = 17° to represent the if-band feature 
of NGC 4501. The amplitude of the deviation from the 
axisymmetric part is taken to be 

*p(^)-^(]^^^ (6) 

which was used for bar models by Sanders (1977) and 
Wada (1994). We solve the equations of motion for test 
gas particles moving in this potential (see appendix 1 for 
the detailed description for damped-orbit model in a spiral 
potential). The resultant orbits are shown in figure 10. 
Oval orbits appear, whose major axes change with the 
radius. As a result, spiral-shaped orbit-crowding regions 
are generated. 

4-2. Comparison with CO Distribution and Velocity 
Field 

Figure 11 shows the density distribution and velocity 
field obtained from orbit calculations. They are projected 
with the inclination and position angle of NGC 4501. The 
density distribution was simply estimated based on the 
assumption that the staying time of a gas cloud at each 
point is inversely proportional to the velocity. The ob- 
served radial profile was applied. After that, we con- 
volved the model images with a Gaussian beam of FWHP 
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Fig. 10. Particle orbits in the fixed spiral potential. An or- 
bit-crowding region is generated as a spiral pattern. Used 
parameters: 8 = 17°, e = 0.03, A = A/k = 0.1, Qb = 0.01, 
a = l. 

65" to reflect the primary-beam attenuation of NMA. Our 
model (figure 10) shows that the clouds slow down near to 
the major axes of the oval orbits, which make the orbit- 
crowding region. The cloud distribution becomes similar 
to the orbit-crowding pattern, and gas clouds moving in 
a spiral potential gather along the spiral pattern. This 
naturally explains the observed molecular arms. The dis- 
tortions along the spiral arms and in the nuclear region of 
the model velocity field well represent that of the observa- 
tions. The velocity field in figure 11 shows that the trends 
in the linc-of-sight velocity of gas clouds change across the 
arms. On the far side of the arms to us, gas is approach- 
ing compared to circular motion. When we look nearer 
and nearer, its line-of-sight velocity approaches that of 
the circular motion. After that, it turns to receding mo- 
tion, and gradually returns to the approaching motion. 
This behavior is well reproduced by this model. 

5. Discussion 

5.1. Central Double Peaks 

Figure 7 shows that the double peaks inhabit the root 
of spiral dust lanes penetrating into the nucleus. The nu- 
clear spiral dust lanes have recently been found in many 
galaxies, and are thought to be the location of shock, and 
trace the sites of angular-momentum transfer, since shocks 
cause a loss of kinetic energy, (e.g. Regan, Mulchaey 1999; 
Martini et al. 2003). The CO "twin peaks" are often found 
in barred galaxies (e.g. Kenney et al. 1992). The "twin 
peaks" are connected to large-scale dust lanes, and regu- 
larly accompany to nuclear star- forming rings. They show 
a gas-to-dynamical mass ratio of as high as ~ 0.3, and 
are thought to be gravitationally unstable enough for star 



formation. They arc often considered to be a the orbit- 
crowding region with an xxjxi orbit population change 
near to the inner Lindblad resonance. On the other hand, 
the gas-to-dynamical mass ratio within r ~ 5" (390 pc) is 
^ 3.5% and no nuclear star-forming ring is seen near to 
the double peaks in NGC 4501 (see figure 5 in Koopmann 
et al. 2001). This is the first time that this kind of "twin 
peaks" was found in a non-barred galaxy, although the 
physical conditions arc significantly different from that of 
the "twin peaks" in barred galaxies. An inner Lindblad 
resonance also exists in the spiral potential, and the ap- 
pearance of double peaks might indicate that they are 
associated with the inner Lindblad resonance, as well as 
those in barred galaxies. However, it cannot be judged 
from our present data only, because the central rotation 
curves are derived from a 5" resolution velocity field, while 
the double peaks are at r ^ 2"4. In fact, the epicyclic fre- 
quency, K = y^{(P^o/dR'^) + Sn^, is sensitive even to small 
errors in the rotation curves. Observations at higher res- 
olution and sensitivity still remain as future work. 

5.2. Central Condensation and Mechanisms of Angular 
Momentum Transfer 

As revealed by Sakamoto et al. (1999a), central CO con- 
densations within r ~ 500pc are prevalent in CO-luminous 
galaxies, and NGC 4501 is one of those with central CO 
excess. In terms of its radial profile, the much shorter e- 
folding radius (4."3 ~ 340 pc; this work) than that of the 
global CO disk (47"0 - 3.8 kpc; Nishiyama et al. 2001) ev- 
idences some distinct origins of the central gas disk from 
that of the global gas disk. Thus, it indicates a possible 
radial inflow of the gas toward the central disk. 

There are two possible mechanisms to transfer angular 
momentum and transport the gas into the central region 
in a global perturbation, like bars or spirals. In the follow- 
ing, we consider these mechanisms in a spiral potential, 
and estimate the effectiveness. 

(1) Galactic Shock Wave: The calculation shows that the 
orbits are crowded along the spiral pattern. It is likely 
that galactic shock occurs in these orbit-crowding regions, 
and that this leads to the loss of angular momentum of 
the gas clouds. We estimated this loss by considering the 
oblique shock of isothermal gas. A detailed description 
is shown in appendix 2. The strength of the shock de- 
pends on the angle between the shock front and the gas 
streamline. With simple approximations, the loss of an- 
gular momentum is ^ sin^ 6 times the initial value after a 
shock. Adopting a pitch angle of 6* = 17° for NGC 4501 
derived from TiT'-band image, the loss is estimated to be 
16% after an orbital rotation. 

(2) Drag of Stellar Spiral Potential: Gas clouds show ellip- 
tical orbits in a nonaxisymmetric potential made by stel- 
lar spiral arms. Therefore, when a cloud is in the leading 
side of the potential valley, gravitational torques from the 
stellar spirals are exerted to drag the cloud 'backwards', 
which result in angular-momentum transfer from the gas 
to the stars. Note that gas clouds will oppositely obtain 
angular momentum when they are on the trailing side. 
Whether a gas cloud loses or obtains angular momentum 
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Fig. 11. Gas distributions (left panels) and velocity fields (right panels) obtained by observations (upper panels) compared to those 
by our models (lower panels). The parameters used for the model arc the same as that used for figure 10. The model disk is inclined 
with the inclination and position angle of NGC 4501. 
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at a certain radius is determined by integration of this 
torque along the orbit of the gas, as shown in equation 
(A23). In a bar potential, the phase of the major axis of 
elliptical orbits always leads the bar potential inside of the 
corotation resonance (CR); therefore, clouds stay longer in 
the leading side of the potential valley than in the trailing 
side. This leads to a loss of the angular momentum and 
inflow of the gas inside CR (Wada 1994). We calculated 
the change of angular momentum, AJz{Ro), for a spiral 
potential of 6' = 17°, and obtained the loss of the angular 
momentum, /S.J z{Rf)) / J ^{Ro) = 0.8%, at most, using the 
parameters presented with figure 10. This value depends 
on spiral perturbation strength, e, and pitch angle, 0. 

The above estimations show that the angular- 
momentum variation due to galactic shock (16%) is 
an order-of-magnitude larger than that of gravitational 
torques (0.8%). Therefore, gas clouds are destined to lose 
angular momentum due to shocks regardless of a gain or 
loss due to the torques from the stellar spiral potential. 

6. Conclusions 

We have observed CO(l-O) emission from the non- 
barred Seyfert 2 galaxy NGC 4501 to study the gas dy- 
namics in the central 5 kpc region. 

1. In the central 5 kpc, NGC 4501 has two molecu- 
lar gas components — spiral arms penetrating into the 
nuclear region, and a central concentration with radius 
5"(- 390 pc), and the gas mass Mgas = 1.3 x 10^ M©. 

2. The spiral arms are observed from the nuclear region 
to the end of our field-of-view (r ^ 40"). The velocity field 
along these arms shows deviations from circular motion. 
That is about ^ 50kms~^ in linc-of-sight velocity on the 
galactic minor axis. 

3. The central concentration is resolved into double 
peaks separated by 4"7(^ 370 pc) in the high-resolution 
map. The double peaks are located on the root of nuclear 
dust spirals found by HST. 

4. The well-ordered non-circular motion reveals that 
the gas is driven by density waves in the central region 
of this galaxy, rather than by stochastic processes. This 
result is consistent with that inferred from near-infrared 
observations. 

5. To understand gas motions in NGC 4501, we studied 
the gas cloud orbits in spiral potential with the damped- 
orbit model. The result well explains the velocity field. 
The gas dynamics in NGC 4501 can be well understood as 
being governed by a spiral potential. The observed density 
distribution can also be explained by a simple assumption 
that the staying time of clouds is inversely proportional 
to its velocity at each point. 

6. Centrally-condensed double peaks have a low star- 
forming efficiency, which may come from gravitational sta- 
bility of gas, owing to low gas-to-dynamical mass ratio. 
It may be possible that they are the results of angular- 
momentum transfer due to shocks in nuclear dust spirals. 

7. Based on our model, we estimated the loss of angular 
momentum in NGC 4501 due to two possible mechanisms 
— galactic shock in orbit-crowding regions and gravita- 



tional torques exerted by the stellar spiral potential. The 
effect of galactic shock is an order-of-magnitude larger 
than that of gravitational torques. Stellar spiral arms are 
responsible for the angular-momentum loss of gas, and 
may lead to gas inflow toward the galactic centers. The 
gas inflow rate depends on the galactic morphology. Our 
spiral model and its consequential gas inflow are natural 
extensions of those discussed for a bar. The spiral-shock 
inflow mechanism, as well as a bar-induced inflow, can 
be a potent instrument to study the relation between the 
galactic morphology and the central gas features, such as 
nuclear condensations and starbursts. 
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ments. We also appreciate Keiichi Wada for his profitable 
comments. S. O. and J. K. are financially supported by 
the Japan Society for the Promotion of Science (JSPS) for 
Young Scientists. 

Appendix 1. Gaseous Orbits in Spiral Potential 

In section 4, we used the damped orbit model to de- 
scribe gas motions in a spiral potential. This model fol- 
lows Binney and Tremaine (1987), which described ana- 
lytic solutions of equations of motion for collisionless par- 
ticles in a weak bar, and obtained stellar orbits under 
an epicyclic approximation. Wada (1994) introduced the 
damping force term considering the collisional nature of 
gas, and obtained gas orbits in a bar potential. Wc modify 
this dampcd-orbit model to introduce a spiral potential, 
instead of a barred potential. We will consider the gas 
clouds moving in a fixed stellar potential. The self grav- 
ity of gas is not taken into account. 

We assume that the pattern of a potential rotates with a 
certain pattern speed, fip. Let {R,ip) be polar coordinates 
in a frame that rotates with the potential, ^{R,tfj). The 
equations of motion of a test particle in this frame are 
written as 



R-Rij^^-^ + 2Ripn„^ 
oR 



^;r, 



R.ij + 2Rij ■ 



1 a$ 



2Rn^ 



(Al) 



(A2) 



Here, we derive analytic solutions of equations of motion 
to the first order. Wc assume that the non-axisymmetric 
part of the potential, $1, is much weaker than the axisym- 
metric part, $oi then 



$(i?,V) = *o(i?) + *i(i?,^), 



(A3) 



where j^i/c&ol ^ 1- We divide R and ip into zeroth- 
and first-order parts: 
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R{t) = Ro + Ri{t); i^{t)=Mt)+Mt)- (A4) 



We then substitute these expressions into equations 
(Al) and (A2). Now, the first-order terms become 



Ri 



dR^ 



Q' 



Ri - 2i?,o^o^'i 



Ro 



a$i 



dR 



-Ro 



(A5) 



i'i = - 



2f7o 



Cf + CI mRn{no ^ np) 
x{{Cifo-C2fi)sm[mino-^p)t + a] 
-{C2fo + Cih)cos[mino-%)t + a]} 



where 



mRlino-np)- 



■ sin[TO(r2o — ^p)t - 



(A14) 



where 



^^ + 2^0-^ = 
Ro 



n{R) = ± 



7?2 



"9^ 



Ho 



1 d$o 
R~dR 



r!o = n{Ro) ; ^0 = («o - ^p)i, 



Ha : 



dR^ 



3f72. 



(A6) 

(A7) 
(A8) 

(A9) 



We now introduce a spiral-shaped perturbation. 
$i(E,V) = $p(i?)cos[m(r!o-^p)i+plog(i?)], (AlO) 

which represents an m-armed logarithmic spiral with p = 
TO cot 0, where 6 is the pitch angle of the spiral arms. Next, 
we add a damping-force term to the equation of motion 
(A5) following Wada (1994). The damping force term 
represents the effects of the collisional nature of gas; thus, 
a plausible approximation can be the form of 2Ai?i, which 
is proportional to the velocity of the radial oscillation, and 
an appropriate value, A. 

Then, the equations of motion become 

i?i + 2Ai?i + Kgi?! =/o cos[m(rio - ^p)^ + a] 

+ fism[m{no-^p)t + a], (All) 



tpi + 2f}o 



Rn 



5$p(i?) 



dR 



cos [m(r2o ^ ^p)t + a] 



Ci 



<Pp{R)-^sm[m{no-np)t + a]. (A12) 

Rn 



We assume the axisymmetric part of the potential, 
^q{R) to be Plummer potential (5), and the strength of 
the non-axisymmetric part, $p(i?), to be equation (6), 
following Sanders (1977) and Wada (1994). When we as- 
sume the gas orbits to be closed, the solutions for these 
equations are 



i?i = 



1 



Cf + Cl 



{{Cifo - C2fi)cos[m{no - np)t + , 



-{C2fo + Cih)sm[m{nn-np)t + a]}, (A13) 



a = p\og{Ro) 



2r2$r 



d$r 



dR R{n-np) 



Rq 



h = $p(i?o) 



da{R) 
dR 



(A15) 
(A16) 
(A17) 



-Ro 



Ci 



''{VLo-VLpf , C2 = 2\m{na-np). (A18) 



Particles trace oval orbits whose position angle of the 
major axes changes continuously due to both the damping 
term and the spiral potential. This situation is different 
from that of stellar orbits in a bar potential, whose major 
axes can only be parallel or perpendicular to the bar. 



Appendix 2 
Variation 



Estimation of Angular-Momentum 

A2.1. Galactic Shocks 



Behaving as a fluid, interstellar gas becomes supersonic 
by a variation of the local gravitational field due to the 
density wave (Fujimoto 1968a; Roberts 1969). This re- 
sults in a galactic shock wave. It is possible that gas loses 
angular momentum due to galactic shocks along the spi- 
ral arms. We estimate the effect of the shock with simple 
approximations of isothermal gas with a Mach number of 
M ^ 1. Although the orbits arc non-circular, a circular- 
orbit approximation is enough for an order-of-magnitude 
estimation. The Mach number, M , is defined as M = v/a 
with velocity v and sound speed a. For interstellar gas, 
because M ~ 10 with v ~ lOOkms"^ and a ^ lOkms"^, 
A# 3> 1 is plausible. 

We now consider the gas bumping into a shock front 
with velocity vi and angle a. In the jump condition of 
isothermal oblique shock, the velocity perpendicular to 
the shock front is reduced by l/A/^. It leaves then there 
with velocity V2 and angle /3, which are written as 



V2 = Vi < 



/sin a 
i\/4 



- cos^ a. 



(3 = tan 



_^ /tana 



(A19) 



(A20) 



Thus the change in the angular momentum, Ji — J2, be- 
comes 
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Ji — J2 = rvi cos 

= rvi 



a) — rv2 cos{6 — /3) 



1 



cos(6' — a) - 



tan^a 



M4 

tana 
172" 

Then, the supersonic condition l/A/^ ^0 and the circular- 
motion condition a ^ 9 yield 



Ji — J2 ^ sin 6* Ji . 
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